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Abstract: 

We use the Inverse Scattering Transform and a diffusion approximation 
limit Theorem to study the stability of soliton components of the solu- 
tion of the nonlinear Schrodinger and Kortcwcg - dc Vries equations under 
random perturbations of the initial conditions: for a wide class of rapidly 
oscillating random perturbations this problem reduces to the study of a 
canonical system of stochastic differential equations which depends only on 
the integrated covariance of the perturbation. We finally study the problem 
when the perturbation is weak, which allows us to analyze the stability of 
solitons quantitatively. 
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1. Introduction 

The aim of the present work is to study the stability of the soliton components 
of solutions of completely integrable systems under rapidly oscillating random 
perturbations of the initial condition. We will consider and compare two impor- 
tant examples of equations widely employed to model nonlinear and dispersive 
effects in wave propagation: the (1-dimensional) nonlinear Schrodinger (NLS) 
equation 

f + iS + »^ = » 

and the Kortewcg - de Vries (KdV) equation 

dU r . TT dU d 3 U n 

The NLS equation models in particular short pulse propagation in single-mode 
optical fibers (then t is a propagation distance and a; is a time) [MN] . The KdV 
equation models shallow water wave propagation [Wh]. 
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Explicit computations are developed in the case of a square (box-like) initial 
condition perturbed with a zero mean, stationary, rapidly oscillating process 
v e (x) := v(x/e 2 ): 



but they can be extended to the case of perturbation of a more general ini- 
tial condition defined by a bounded, compactly supported function q(x). Our 
approach to both examples relies on the results of the Inverse Scattering Trans- 
form (1ST), a powerful tool to study solutions of completely intcgrable nonlinear 
equations, see [APT]. In this framework, the problem is transformed into a linear 
system of differential equations where the initial condition enters as a potential 
and soliton components correspond to eigenvalues. Indeed, any initial condition 
may generate soliton components, that are solitary waves that propagate over 
arbitrarily large distances with constant velocity and constant profile, and ra- 
diation components, that decay in amplitude as a power law. The identification 
of the soliton components therefore characterizes the long-time behavior of the 
solution of the PDE. 

We will show in Section 2 that for rapidly oscillating processes (small values 
of e) the limit system governing the stability of the soliton components reads as 
a set of stochastic differential equations (SDEs) and it is formally equivalent to 
the system where the initial condition contains a white-noise perturbation: 



where a is the integrate covariance of the process v. This shows that to study the 
soliton components in the limit of rapid oscillations, the only required parameter 
of the statistics of v is its integrated covariance. Notice that we cannot directly 
use a white noise to perturb the initial condition, as the 1ST requires some 
integrability conditions on the initial condition (for example, Uq £ L 1 ), which 
are not satisfied by a white noise. 

Our results state that solitons are stable under perturbations of the initial 
condition for both examples studied. However, a few interesting differences will 
be pointed out; in particular, thresholding effects for the creation of solitons are 
present in the NLS case and absent for KdV. 

We also provide an easy way to compute the first order corrections to the 
parameters characterizing the soliton components of the solution. 

In the 1ST framework, a direct scattering problem (known as the Zakharov- 
Shabat spectral problem - ZSSP) associated to the NLS equation is introduced 




(3) 




(4) 



(5) 



where ip n are the components of a vector eigenfunction ^ and £ is the spectral 
parameter. For any localized initial profile Uq the continuous spectrum consists 
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of the real axis. For ( G K the solutions of (5) with the boundary conditions 

* ~ (J) e ~ KX * ~ (l) ' x^-oo (6) 

define two sets ('J', 'J') and (<&, $) of linearly independent solutions. Therefore, 
these functions arc related through the system 

*\ = (KO a(0\ (* 

*y \ho ho) \<£ 

where a, b are called Jost coefficients. The functions ^ and $ are called Jost 
functions. If Uq E the function a(£) is analytic in the upper half of the 

complex plane, where it can only have a countable number of simple zeros, see 
[APT, Lemma 2.1]. These zeros are the eigenvalues of the discrete spectrum. 
Each discrete eigenvalue £ = £ + if], r\ > 0, corresponds to a soliton component 
of the solution. A pure soliton solution of the NLS equation has the form 



exp 

U(t,x) = 2iri 

cosh ( 2i](x + 4£i) 
up to a shift and a phase. 

Similarly, one can link the existence of soliton components of solutions of 
the KdV equation to the spectral properties of an associated equation, the first 
equation of the Lax pair: 

g + (c/„ + C 2 V = o, (7) 

where <p : M. — > K.. We need to assume that 

|/:M^M f (1 + \xf)\U(x)\dx < oo} (8) 



Uo e P„ 

for fi = 1 or /i = 2, see [AC, chapter 2]. Consider the continuous part of the 
spectrum of equation (7), which is again the real axis. For (el, there are 
two convenient complete sets of bounded functions associated to equation (7), 
defined by their asymptotic behavior: 

4>{x, C) ~ e~ iCx , 4>{x, C) ~ e Kx for x -> -oo 

t/)(x, C) ~ e lCx , tp(x, C) ~ e~ lCx forx^+oo. 

It follows from the above definitions that 



<Kx,C) = <Kx,-£), il>(x,0=i>(x,-C) 
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and 

^,0 = 0(0^,0 + 6(0^,0 

0(aj,O = -o(O^»O + KO^ ) O ■ 

The function a is analytical in the upper half of the complex £-plane, where it 
has only a finite number of simple zeros located on the imaginary axis f = irj, 
see [AC, Lemma 2.2.2]. These zeros are the eigenvalues of the discrete spectrum, 
and they correspond to the soliton component of the solution. A pure soliton 
solution is given by 

U(t, x) = 2rj 2 sech 2 (77(31 - x(t))) , 
where x{t) = xq + Ar) 2 t is the center of the soliton. 



2. Limit of rapidly oscillating processes 

This section contains a rigorous justification of the use of the 1ST. As remarked 
in the introduction, to be able to apply the 1ST, the initial condition Uq needs 
to satisfy some integrability condition, L 1 for NLS and (8) for KdV. These 
hypotheses are satisfied by initial conditions of the form (3) if v is bounded. 
Our objective is to show that the 1ST applied to these random initial conditions 
gives a problem that reads as a canonical system of SDEs in the limit e — > 0. 
We make the following assumptions on the process v. 

Hypothesis 2.1. Let v{x) be a real, homogeneous, ergodic, centered, bounded, 
Markov stochastic process, with finite integrated covariance J Q °° E[z/(0)i/(x)]dx = 
a < 00 and with generator C u satisfying the Freedholm alternative. 

Set 

Uq ■= (<z + %(x/e 2 ))l M (x) 
and remark that for 16 [0, R] 

\ U^(x')dx' ^> / U (x')dx' = qx + V2^aW x 
Jo Jo 

in distribution, see [FGPS]. For every e > 0, we apply the 1ST to the NLS and 
KdV equations with the initial condition Uq. 

We consider the ZSSP associated to the NLS equation: our goal is to identify 
the points of the upper half of the complex plane, ( € C + , for which there exists 
a solution VP of the first order system (5) for x € [0, R], satisfying the boundary 
conditions 

*(0)=(T) and ViOR) = 
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derived from the exponentially decaying conditions (6). These particular values 
of ( are the discrete eigenvalues of the ZSSP and correspond to the soliton 
components. The strategy employed is to consider the flow ^(x, £), x <G [0, R], 
C € C + , solution of (5) with initial condition 

*(0) = (J) (9) 

and look for the values of C for which the final condition is satisfied. 

For a fixed value of £, we consider the solution ^ £ of the ZSSP obtained from 
the 1ST 

m = - iCri+i us(x)r2 
^i = t (u^x)yr 1 + Kr 2 



with initial condition 



* e (0) = 



Now, [FGPS, Theorem 6.1] states that the process VP 6 converges in distribution 
in C([0, R]; C 2 ) to the process \P solution of 

{dipi = [(— i( — aa 2 ) ?Ai + iqip2\dx + iy/2a afodWx 
dip2 = [iqipi + («C ~ aa 2 ) ifa] dx + i\f2aail>idW x 
with initial condition (9), which can be rewritten in Stratonovich form as 

d* = i (~G f] $dx + iV2a~ a ( J * ) * o dW x . (12) 



q CJ ^ 

For NLS we can also consider perturbations produced by a complex process: 
let v\,V2 be two independent copies of the process v and set v := V\ + iv%. 
One can define Uq using v instead of v; proceeding as above, from the 1ST one 
obtains again the system (10), and from [FGPS, Theorem 6.1] one gets that in 
this case the limit process is the solution of 

d^=i(~^ ^¥da + *V2acrfj Jj^odW^-V^ao-f^j J J * o dWj 2) , 

. a . 3) 

where the arc two independent Wiener processes, with the same initial 

condition (9). 

We apply the same strategy to the KdV equation: the goal is to obtain the 
values of £ € C + for which there exists a solution <p of 

^ + W + CV = (14) 
with the boundary conditions 

¥>(0) = 1, ^(0) = -iC, <p x (R)-iC<p(R) = 0. 
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These conditions correspond to imposing exponential decay of the solution at 
infinity. 

Setting $ := (ip, tp x ) T this equation can be transformed into 

d$ £ = (_°_ c2 J) ^dx , $(o) = , MR) - = o • 

We consider the flow & £ (x, C), x E [0, R], £ € C + , defined by the above equation 
with only the initial condition, and look for the values of £ s.t. the final condition 
is satisfied. Again by [FGPS, Theorem 6.1], $ £ converges in distribution to the 
solution of 



\_ q °_ C 2 jWx + Vwr ° Q ) 



which, in terms of the function tp, can be rewritten as 

d<p x =-(q + C 2 ) <pdx + V2aaipdW x . (15) 
The initial condition is 

$(0) = > or equivalently (p(0) = 1 , ^(0) = -< • (16) 

Remark that in the last two differential equations above the Stratonovich and 
ltd stochastic integrals coincide. 

The convergence obtained above is only for a (finite number of) fixed £ and a , 
but we will need a convergence in C°([0, R]; C 1 (M 3 )) to be able to differentiate 
the limit process with respect to the parameters. This is the main result of 
this section and it is provided by the following Theorem. We will focus on the 
problem of finding the values of £ for which the limit flows W and $ match the 
final conditions in Sections 3 and 4. 

Theorem 2.2. Assume Hypothesis 2.1. Let "J 6 := {4't^2) T be the solution of 
(10) with initial condition (9) and the solution of (12) with the same initial 
condition. Let also ip E be the solution of (14) with initial condition (16) and 
tp be the solution of (15) with the same initial condition. Considering these 
as functions of the space variable x and the parameters £, n, a, we have in the 
limit of e -> that $> £ {x,£,<q,o-) -> y(x,£,r),a) in C°([0,i2];C 1 (R 3 ;C 2 )) and 
cp s (x, e, f], <r) tp(x, f , J], a) m C°([0, R}; C 1 (R 3 ; C)) . 

To prove this theorem we need the following standard tightness criteria, see 
[Me84]. 

Lemma 2.3. Let (E, d) be a metric space, and X E a process with paths in 
2?([0, R\; E) . If for every x in a dense subset of[0,R] the family (X e (x)) se ^ Q j, 
is tight in E and X £ satisfy the Aldous property: 
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A : For any n > 0, A > 0, there exists S > s.t. 

limsup sup sup F(||X £ (t + 9) - X £ (r)|| > A) < k , 

£->0 T<RO<e<5/\(R,-T) 

where t is a stopping time; 
then the family (^ e ) £g ( ji * s W^/ii in 2?([0, i?]; i?) . 

Lemma 2.4. Let H be an Hilbert space and H n be an increasing sequence of 
finite-dimensional subspaces ofH s.t., for any h £ H, lim.n_j.oo tt-h„ h = h. Let 
X s be an 'H-valued process. Then, X e is tight iff for any k > and A > 0, there 
exist p K and a subspace "H k ,a s.t. 

sup P(||J_ e || > p K ) < k and sup P(d n (X e , U K ,\) > A) < n . (17) 

eG(0,l] £6(0,1] 

Proof of Theorem 2.2. Since Propositions 3.2 and 4.2 ensure that the limit equa- 
tions for and $ have a unique solution which is C°([0, i?];C 1 (K 3 )), it suffice to 
prove convergence in the space of CadLag processes 
P([0, i?];C 1 (R 3 )), which is in turn provided by Lemma 2.3 if we take E to 
be the Hilbert space H := W 3 ' 2 (G), which is imbedded in C 1 (G). Here, G is 
an open, bounded subset of Bl 3 , the space of parameters. For simplicity we take 
G = (-N, N) 3 for some real positive constant N; a justification of the fact that 
it is not restrictive to assume that the set of parameters G is bounded is given 
below in the proof of Proposition 3.2, where the convergences we are proving 
here will be used. 

We start by using Lemma 2.4 to obtain the necessary tightness of the se- 
quences { v I /£ } and {<_ £ }. Let d% be the distance induced on H by the scalar 
product and let X s be the solution of the approximated system (X e is equal to 
either \_ e or Q £ ). The space H n is constructed as follows. Divide G into cubes 
with sides of length 1/n and add one extra layer of cubes around it: Ajj^ := 
[i/n, (i+l)/n)x[j/n, (j+l)/n)x[k/n, (fc+l)/n), for i,j,k = -(nN+1), ...,nN. 
Define the pieccwise (on every cube) polynomials of fourth degree as 

Nn 4 

M*) : = E E — { < a Sl i x - y^) {m) Hm^M . (is) 

i,j,k=-(Nn+l) m=0 

where yij.k is the center of the cube Ai t j t k, o!f"j \ are families of m-dimensional 
tensors and the brackets denote the relative tensor products (so that, for ex- 
ample, (a^jf. \x — Di,j,k)^ denotes the product between the four-dimensional 
tensor afj k and four copies of the vector (x — yi,j,k) )• With this definitions h is 
a function defined on [— N— -\, N+ i) 3 , but its restriction to G does not belong- 
to H in general, since it may not even be continuous. Let T be a real, nonnega- 
tive, smooth function, with compact support contained in [—1/2, 1/2] 3 and such 
that /r_]/2 i/2] 3 r( x ) ^ x = 1- Setting r n (a;) := n 3 T(nx) we can finally define 
H n as the finite dimensional space of functions of the form h(x) := (h*T n )(x). 
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Remark that h has been defined on a set larger than G, so that the convolution 
product is well defined for x £ G. 

Since X £ is the solution of a linear differential equation (if q is constant it 
even has constant coefficients in x) with coefficients smooth in the parameters 
fi = (£,??, er), from the explicit formula for the solution we get that X 6 (x,fi) is 
smooth in the parameters . We will soon use its derivatives in the parameters: 
the vector of X e and its first derivatives in /i still satisfy a linear system of 
ODEs whose coefficients depend linearly on the parameters and on the process 
v(x/e 2 ), and the same result holds adding higher order derivatives. 

The key step to show that {X £ (x, /J*)} &K is tight in H = W 3 ' 2 (G) for every 
x £ [0, R] is the proof of the bound 



limsup E ||X e (a;,-)|| 62 . . < oo (19) 



uniformly in a; £ [0, R]. Indeed, the first part of condition (17) immediately 
follows, and for the second part we only need to use Sobolev's imbedding 
W 6 ' 2 <-» C 4 (R 3 ) and the following lemma. Corollary 2.6 will finally provide 
the last hypothesis of Lemma 2.4. 

Lemma 2.5. For every g £ C 4 (G), there exists a g n £ TL n s.t. 

\\g-9n\\ H < c -h\\ C 4 {G) ■ 

Proof. Step 1: (Construction of g n ). For k — —nN, ... , nN— 1 (which means 
that y hj . k e G) set av?^ := D m g(y l .^ k ). For i,j,k such that y iijtk £ G set 

a i"j\ := D m g(y'), where y' is the nearest cube center; notice that the distance 
of these two points is at most the diameter of the cubes, which we call 25 := 
2 v / 3n~ 1 - With the afj k thus defined we construct the piecewise polynomial 
function lj n as in (18): on the cubes the centers of which are not in G, this 
function is just a copy of the function defined on the nearest cube with center 
in G. Finally, we have g n := T n *g n . 

Step 2: (Estimates). For every multiindcx a £ N 3 such that |o|i := a\ + a,2 + 
«3 < 3, we need to estimate 



d a {Y n *g n ){x)-d a g{x) 



2 

da; 



To clarify the procedure to obtain an estimate for the above term we first give 
explicit computations for the case a = e\ = (1, 0, 0). Recall that, by definition, 

J2 [ r n (x-y)dy= [ T n (x-y)dy=l 
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For a = (1, 0, 0), using twice the inequality (a + b) 2 < 2(a 2 + b 2 ), we have 



di{T n *g n ){x)-d W {x)^= \J2 f -d yi T n (x-y)g n (y)dy - d l9 (a 

i,j,k JAi.i.l* 

< 2{| J2 [ r "( x - y)d yi 9n(y)dy - d ig (x) 2 



/ rn ( x -v) 9n(vt, 2/2, 2/3) -gn(Vi ,2/2,2/3) dy 2 dy 3 



< 



4 {|E / r n (x-y)\d yi g n (y)~d 1 g n (x)\dy 



+ 



dig n (x) - dig(x) 



V / T n (x-y) g n {yt,V2,yz) - 9n (vi ,2/2,2/3) 



dy2 dyj 



4{^ 1 + 5' 2 + 5 3 } 



(20) 



where diAij ^ denotes the faces of the cubes orthogonal to the direction ei := 
(1,0,0). Notice that the number of non zero terms in the sums over i,j, k of this 
proof is limited to 8 because the support of T n can intersect at most 8 cubes. 
The term Si can be bounded by the square of 

ll r "0- ^h^A^^Wdigni-) - dig n (x)\\ L ~ { B^(x)) , 

where the second term is regarded as a function of y (x is fixed) and the L°°~ 
norm is taken on the ball B±_ (x), which is the support of r n (a; — y). The first 
term above is 1, and to estimate the second term the worst case is when y does 
not belong to the same cube as x: let us say that y € A^ 1 ' and x € A^ 2 \ where 
i/ 1 ' and y^ are the centers of the cubes A^> and A^ respectively. We have 
therefore the bound 



\\dig n (-)-dig n (x)\\ L oo iByi2n)(x)) 

< WdM) ~ d yi g n (y {1) )h~ + \d l9 (y {1) ) ~ d l9 (yW)\ 

+ |3 yi 3„(y (2) ) ~ d yi g n {x)\ 
<4S\\D 2 g\\ L ~ {G) . 

This provides the bound for the term Si. Similarly, for S2 we have the bound 

\dig n (x) - d l9 (x)\ < \d{g n (x) - d l9 (y^)\ + \d ig (y^) - d l9 (x)\ 
< 25\\D 2 g\\ L ^ [G) . 

We still need to estimate S3, which contains the boundary terms deriving from 
the discontinuities of g n (and, in the general case, of its derivatives) . This term 
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requires more careful estimates. With Cr '■= sup x T(x) and using the fact that 
l|r™(a; - ^WmdiAij,*) < nC r , we have 



L°°{d 1 A,. ]tk ) 



< nC T Y \ \\Sn(vt 

i,j,k 



■) -g{vu-)\ 



L~(aiA Sj3ifc ) 



+ ||s(W,')-Sn(l/i,0|L« (9lAjJifc) 
< 8nC r 6 4 \\D 4 g\\ L ~ {G} = C6 3 \\D 4 g\\ L ~ {G) . 

The sums over i,j, k above are meant as sums over the faces d\Ai.j^ intersecting 
the support of T n (x — •), which are at most 4. Collecting all these results, we 
have the uniform bound 

9r(r" - d l9 (x)\ 2 < C6 2 (\\D 2 g\\ 2 L ^ {G) + S 4 \\D 4 g\\ 2 L ^ {G} ) . 

We proceed in the same way with higher order derivatives to get, for a generic 
derivative a of order < \a\ < 3, the estimate 

d a (T n *g n )(x)-d a g(x) 2 <C\\J2 I T n {x - y)\d a g n {y) - d a g n {x)\dy 



+ \d a g n (x)-d a g(x) 
<c(\\D^ +1 g\\l x{G) 6 2 + S%) . 

For a = 0, S§ = 0, but in the general case the estimate of the term S§ is 
a little bit more delicate, since one gets more boundary terms. In particular, 
when integrating by parts, derivatives along different directions result in terms 
containing discontinuities of g n and its derivatives along the faces (that we 
denote for brevity dA), edges (denoted d 2 A) or vertices (denoted d 3 A) of the 
cubes, while multiple derivatives along the same direction result in derivatives 
of r™ appearing. For example, for a = (1, 1, 1) we get three kinds of terms: 



L 



T n (x-y) Ag n (y)dy , 



d 2 Ai 



T n (x-y) Adg n {y)dy 



[ T n (x-y)Ad 2 g n (y)dy , 

JdA itj<k 



where A denotes the jump of the function. For a = (3, 0, 0) we also have terms 
like 

/ d 2 r n (x-y)Ag n (y)dy , f diT n (x - y) Ad{g n (y)dy , 
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and in the general case we find also terms like 

dT n (x-y) Ag n (y)dy. 



J 

Jd 2 Ai 



However, we can bound all this terms in the same way. We have that, for m G N 
and 6, c e N 3 , 

/ d b T n (x)dx < C r n m+ l fc l = C(5- (m+|h|) , 

\Ad c g n (y)\ = \d c g n (y + ) ~ d c g n (y-)\ < 2\\D 4 g\\ L ~ (G) 5 4 ~^ . 

Therefore 

SI < C\\D 4 g\\l^ {G) 6 2 ^-^ . 

Remark that for all the terms composing S$ we always have m+|6| + |c| = \a\ < 3. 
Summing up, we have obtained the bound 

\\9-9n\\i=y2 I \d a (T n *g n )(x)-d a g(x) 2 dx 

a J G l 

<C^ (\\D^ +1 g\\U {G) S* + WD^i^S 2 ^) 

a 

<C^\\g\\ 2 c HG) . 

The Lemma is proved. □ 
Corollary 2.6. For any h <E H, lim„_ > . o0 n<n n h = h. 

Proof. Fix any e > 0. By density, there exist a h e £ C A (G) s.t. \\h — h e \\u < e/2. 
Also, by the continuity of the projection, \\7Tn n h — 7Tn n he\\n < \\h— ft. e ||-K < e/2. 
Since \\iru n h e — h e \\-u < \\h £ , n — h e \\u, by the above lemma we get 

\\n-H n h- h\\ H < \\iru n (h ~ h e )\\ n + \\^u n h £ - h e \\ n + \\h s - h\\ n 

<e + C-\\h s \\ CH K). 
n 

Therefore 

limsup ||7r-H n /i — < e , 

n— >oo 

and since e is arbitrary, the Corollary is proved. □ 

We return to the proof of Theorem 2.2. To show that in our case the hy- 
pothesis of Lemma 2.4 are met, we have to show (19). Define Y £ as the vector 
process of X E and all its derivatives in the parameters fj, = (£, rj, a) up to order 
6. As remarked above, this process is the solution of a linear system of ODEs 
with coefficients (the matrices Mi and Mjj) linear in the parameters: 

— Y e = M x Y e + -v{x/e 2 )M 2 Y e . 
dx e 
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Since G is bounded, we only need to check that the second moment of Y £ {x, ji) 
is uniformly bounded with respect to e £ (0, 1], /i £ G and x £ [0, R]. Actually, 
we aim at a stronger result, which we will need later. We are going to show 
that (recall that Yo = Yq is deterministic since it is defined by the equation for 
x < 0, which is deterministic, and the boundary condition at x — > — oo) 



E 



sup |F e (a;)| 
xe[o,R\ 



< C R 1 



r(i + \Y(0)\ 2 ) 



< oo . 



(21) 



Following [FGPS, section 6.3.5], we show this bound with the perturbed function 
method. Let C £ be the infinitesimal generator of the process Y £ and C the 
infinitesimal generator of the limit process Y := lim e ^o Y £ . Let m £ N be such 
that Y e (x) £ C m and let K be a compact subset of R containing the image of 
the bounded process v. For every y £ C m and z £ K, C e has the form 

C e g{ Vl z) = ^jC„g(y, z) + ~ (M 2 y) T V y g(y, z) + (M iy ) T W y g{y, z) . 

Let / be the identity function on C m and f e (y, z) = y + efi(y, z) be the associ- 
ated perturbed function, which is solution of the Poisson equation L v f\{y, z) = 
— z(M 2 y) V y f(y). In this equation, y plays the role of a frozen parameter, so 
that fi has linear growth in y, uniformly in z, and the same holds for 

£ e f%y,z) = z(M 2 y) T V y f 1 (y,z) + e(AI 1 y) T V y f 1 (y,z) . 

Since 



Y £ (x) = Y-{0)-e[h{Y%x)^{x)) - /i(y £ (0), v*(0)) 

+ [ C e f s {Y e {x'),u e {x'))Ax' + M e x , 
Jo 

where MJ is a vector valued martingale, we get the bound 



sup \Y £ (x)\ < \Y £ {0)\ +eC 1 + sup \Y £ (x)\ 

c£[0,R] L xe[Q,R] 



f R 

C 1+ sup \Y e {x')\dx + C sup |M|| 
x'e[o,x] xe[o,R] 



For e < 27^, applying Gronwall's inequality and renaming constants we get 



sup \Y e (x)\ < C R (l + \Y £ (0) \ + sup \M, 

cfn pi xe[0,fl] 



xe[o,R] 



(22) 



The quadratic variation of the martingale is given by 



(M%= / g e (Y s (x'),v E (x'))dx' 
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where 



g E (y,z) = {C e f a -2f e C*f e )(y,z) 
= {C v fl-2hC v h){y,z) 



2ez 
2e 2 



{M 2 yf h{y,z) - {M 2 y) T ((V y f 1 ) T f 1 )(y,z) 
{M l yff 1 (y,z)-(M 1 y) T ((V y f l ) T f 1 )(y,z) 



has quadratic growth in y uniformly in z € K. Therefore, by Doob's inequality, 



E 



sup | Ml 

xG[0,R] 



< CI 



< c 



1 + E y £ (a;) ]dx . 



Substituting into the expected value of the square of (22) and using again 
Gronwall's inequality, we get (21), which gives (19). We can therefore apply 
Lemma 2.4 and obtain the tightness of (A e (x, /")) eif hi % = W 3 ' 2 (G) for ev- 
ery x € [0, R]; thanks to Lemma 2.3, the tightness in 2?([0, R]; E) follows if we 
show that the Aldous property [A] holds. Since G is bounded, we can prove the 
Aldous property showing that 



lim limsup sup sup sup E 

£^0 fi£K t<R 0<6<6 



\Y%r + e,fi)-Y £ (T,fi)\ 



, 



(23) 



where Y £ is the vector process having as components X s and its derivatives in [i 
up to the third order only. We prove the above limit using again the perturbed 
test function method. With the notations introduced above, we have 



\Y e {r- 



Y e (r)\ <C 



M £ T+e - M* 



C 



T + 6 



£ s r(Y s (x),v s (x)) 



Ce(l+ sup \Y £ (x)\ 



< c 



2 



M e T+ e - M: 



C 



dx 



dx 



C 



C c f B {y e ,if)-Cf<Y e ) 
f + \Cf(Y e )\ 2 dx + Ce(l+ sup \Y e {x)\ 2 ) 



Wc have that 



E 



\M. 



Ml 



= E 



= E 



T + 



d(M<% 



{M s T+e ) - (M r e ) 

Since \£\r(y, z) - Cf(y)\ < eC(l + \y\) and \Cf(y)\ < C\y\, for 9 < S 
e[|F £ (t + 9) - Y e (t)\ 2 ] < C R (S + e)(i+e\ sup |y £ (.T)| 2 ]) . 

L - 1 ^ ^-xe[Q.R] ^' 

The right-hand side is independent of r and we can use estimate (21) to bound 
it uniformly in e and fi. Therefore, (23) follows and the proof of Theorem 2.2 is 
completed. □ 
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Remark 2.7. As seen in [EK, Theorem 7.4.1], the conditions imposed on the 
driving process v(x) can be relaxed, assuming that it is just a mixing process. 



3. Stability of NLS solitons 

This section is devoted to the study of our first example, the NLS equation. 
We focus on the soliton components of the solution. In the previous section 
we have obtained the limit equation (12) and we have seen that every soliton 
component (soliton, in short) is identified by a complex number £ = £ + in 
s.t. the flow ^(x, £) solution of (12) with initial condition (9) satisfy also a 
given final condition. The real and imaginary parts of £ define the velocity and 
amplitude of the soliton, respectively. We start by reporting, in Paragraph 3.1, 
some classical results on the background deterministic solution. Then we analyze 
how this solution is modified by the introduction of a real, in Paragraph 3.2, 
or complex, in Paragraph 3.3, small-amplitude white noise perturbation of the 
initial condition. The main results are contained in Propositions 3.1 and 3.6. We 
deal with the limit cases of "quiescent" solitons in Corollary 3.4 and Remark 
3.8. 

For simplicity of exposition, in the present and following sections we will 
choose the value of the integrated covariance of the process v to be a = 1/2. 

3.1. Deterministic background solution 

Let Uq(x) = q lro.ra (x). Burzlaff proved in [Bu88] that in this case the number 
of solitons generated is the integer part of 1/2 + qR/ir (see also the relevant 
discussion and generalization of [Ki89]). They remark that physical intuition 
suggests that the first soliton created when increasing R corresponds to £ = 
(this is a single soliton with zero amplitude and velocity, the quiescent soliton); 
this "soliton" is created for qR = ir/2. For values of qR just over this critical 
threshold the created soliton has zero velocity and nonzero amplitude 2n which 
can be computed explicitly solving (5) for pure imaginary values of £. 

In the first part of this paragraph we report some computations relative to this 
case, as the results and explicit formulas will be used below. We then conclude 
the paragraph providing the sketch of an analytical proof of the claimed fact 
that generated solitons correspond to purely imaginary values of £. 

From the decaying condition at — oo one obtains the initial condition 

*(0) = (J) e 
The system (5) for x <G [0,R] reads 

^ = iqip2 - Ki>i 



(24) 
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and $ = {ipi, 1P2) is a solution of the initial value problem for £ 7^ if 

V^+c 2 



= 7^ i a 2 sin g2 + ^ x ) + c °s(vV + C 2 a?) ; 



Mx) = « r^— s sin (vV + C 2 x) ■ (25) 
V 9 + C 2 

To be a soliton solution, vff needs to satisfy also the decaying condition at +00, 
which is to say tpi(R) = 0. The condition can be rewritten for £ ^ as 

/ = tantv/^TC 2 ^) + j VV + C 2 = _ (26) 

Since a(£) = ipi{R, C)e ll ' R , the function / is linked to the first Jost coefficient a 
by the relation 

e -iCR y/gJTC 



f = ia(0 



C COS ( V / q 2 +C 2 R) ' 



from which we see that the zeros of / coincide with those of a. 

Note that for £ = iq the solution of (3.1) satisfying the initial conditions is 



1 + x 

ix 



which cannot satisfy the final conditions, so that no soliton can be created for 
this particular value of £. 

Here, plotting on the (£, rj)— plane the level lines of the real and imaginary 
parts of / at height zero, one can already see that the first soliton component 
of the solution corresponds to an imaginary value of £. This fact can be proved 
in an analytic way as follows. 

Recall that the zeros of / coincide with those of a and observe that the 
function a(£, rj, R) is analytic in the domain K x (0, 00) x (0, 00) and continuous 
in M. x [0, 00) x (0, 00). We use the argument principle to study how the number 
of zeros in the upper half of the complex plane evolves with increasing R. For 
any fixed R we proceed as in [DP08], taking a loop C in the complex £-plane 
composed of the (lower) real axis and the infinite semi-arc in the upper half 
plane. Then, the number of of zeros is given by 



2tt Jc a dC, 

Since a = 1 + 0(l/() for \(\ ^> 1, the integral over the upper part of the loop is 
zero. Changing variables a(£) = p(() exp(ia(C)), after some computations one 
obtains that unless there is a zero on the real axis, also the integral on the lower 
part of the loop is zero. Therefore, the number of zeros changes for a given R 
only if a(£, 0, R) — admits a solution. But zeros of a and / coincide, and since 
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in equation (26) for real values of ( = £ 7^ the tangent is real and the second 
term is purely imaginary and non zero, solutions of /(£,0, R) = can only be 
found at £ = 0. 

Explicit computations easily show that C = corresponds to a soliton solution 
only for R = 2n ^ q l 7T, n£N. Computing explicitly the derivative of a(() at £ = 
we get 



9 c a(C) 



CR 



d c a(() 



i— sm(qR) 



+ iR- 



sin 



r 



iRcos(qR) 



so that for Rq = 2 "+ 1 ir the derivative is equal to ^fi/q and is never zero. 
Therefore, new solitons are generated one at a time and they are immediately 
pushed (as R increases) towards the interior of the domain. [Ka79] showed that 
if a(£) = then a'(£) ^ 0; from this fact it follows that zeros in the interior of 
the domain are always simple. Considering the complex conjugate ^/*, which is a 
solution whenever $ is, one obtains that zeros not laying on the imaginary axis 
always come in pairs ±£ + in. But since zeros move continuously (as R grows) in 
the upper complex plane, cannot coalesce and cannot leave the imaginary axis 
(£ = 0) unless they form a pair, we get that they must remain on the imaginary 
axis. 



3.2. Small-intensity real white noise 

In this paragraph we consider the example of an initial condition composed 
of a square function perturbed with a small real white noise. First, we use a 
perturbative approach to study the effects of the perturbation on the first soliton 
(Proposition 3.1); the effects on subsequent "true" solitons are the same. Then, 
in the last part of this paragraph, we study the effect of this perturbation on 
"quiescent" solitons (Corollary 3.4). 

We have here Uq(x) = (q + aW x ) 1[ 0i _r](x). The initial condition is (24) and 
the system (5) for x € [0, R] reads 

{dijji = i{qip2 - C V'l) da; + ia ^2 dW x 
(27) 
d^ 2 = i(qipi + C ^2) dx + iaipx o dW x 

Proposition 3.1. For qR > -| and in the limit of a small real white noise 
type stochastic perturbation of the initial condition, the amplitude of the soliton 
component of the solution is perturbed at first order by a small, zero-mean, 
Gaussian random variable, which is given by 

q sm(c R) W R + f* 222f sin (c (R - y)) sm(c y)dW y 

! ( 28 ) 



% + Rr] sin (coR) - R^ cos (c R) 
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where cq := y/q 2 — T]q, and r/o is the amplitude of the soliton of the unperturbed 
system. 

The velocity of the soliton remains unchanged. 

The following proposition provides results that we will need in the following 
on the uniqueness and regularity of solutions. 

Proposition 3.2. The stochastic differential equation (27) defines a stochastic 
flow ^x ,a ^ = {%j}i,ip2) T of C 1 -difjeomorphisms, which is C 1 also in the parame- 
ters £, rj, a . 

Proof of Proposition 3.2. Write the SDE in Ito and vector form: 

d*=*(- c ;^ 2 c+ q j)**°+i*$ (29) 

The coefficients of the SDE are independent of x and Lipschitz continuous in 
for every £ and a. Therefore, the existence of a unique solution to the SDE, 
which defines a stochastic flow of homeomorphisms (&), is a classical fact 

(see for example [Ku84]). Following the notation of [Ku], we define the local 
characteristic of the SDE as (a, 6, x), where 



1 



b(C, <r,x) 



-C q 
q C 



Fix any n £ N, define the set G n := {(£, 77, <r)| |£| < n, < r\ < n, < a < 
n} and consider the SDE only with parameters in G n . Then, both a and b 
are uniformly bounded and, together with their first derivatives, are Lipschitz 
continuous in the parameters. This means that the coefficients satisfy condition 
(A5)i )0 of [Ku, Chapter 4.6]. It follows from [Ku, Theorem 4.6.4] that 
is C 1 in the parameters almost surely on G n . Let fl n be the set of w € f2 such that 
g C 1 (G n ): it is a set of full measure. Since n is arbitrary, ^^ ,a >(x) 
is actually C 1 (R x R + x K + J for every u> £ n n f2„, which is still a set of full 
measure. This proves the last statement of the proposition. 

Since [Ku, Theorem 4.6.5] states that ^^'"^(x) is actually a stochastic flow 
of C 1 -diffeomorphisms, the proof is completed. □ 

Proof of Proposition 3.1. Looking at the flow at point R we can define a complex- 
valued function of ^(R) as F(£, r), a) := ip^'^ (R). We look for the set of values 
of (£, rj 1 a) corresponding to zeros of the function F: they are the parameters 
(( = £ + in) defining the soliton components of the solution of the problem 
perturbed with a noise of amplitude a. We claim that a small stochastic per- 
turbation has only the effect of a small variation in the value of £ = £ + in with 
respect to the value Co = i-Vo °f the corresponding soliton in the deterministic 
case. We will prove this using the implicit function Theorem: F(£, 77, a) has a 
unique zero in some open set containing the point (0, rj , 0). The following lemma 



E. Fedrizzi/ Random perturbations of solitons 



18 



ensures that the most important hypothesis of the implicit function Theorem is 
satisfied. 

Lemma 3.3. Let F(£, 77, a) be the function defined above. Then, whenever £0 = 
ii]Q is the value corresponding to a soliton component of the solution of the 
deterministic problem, the determinant of the Jacobian matrix 



J 



d v 3(F) J 



at point (0, 770,0) is not zero. 

Proof of Lemma 3.3. For a = the system (27) becomes deterministic and the 
solution is given by (25). Then, setting c 



so that 



r g2 



iR 



sin (cR) 



c 2 

R— cos 

c z 



cR) 



1 %F(0, 770,0) = d v F(0, 770,0) 



We are left to verify that 
7^ det J (0,77 ,0) = 



9 e n(F)d n 3(F) - d n 3£(F)% 3(F)J (0, r? , 0) 
R(0 £ F)Z(d v F) - n(d v F)$(d 6 F)] (0, 770, 0) 
(5R(^F)) 2 + (»(^F)) 2 ] (0,770,0) 



or equivalently 

%F (0,r7o,0) 



eg c 



sin (c i?) + cos (c R) ^ 



(30) 



where Co = \/ q 2 — rfe . Observe that condition (26) implies that 770 < q; cq is 
therefore real. Indeed, for 770 > q, Co would be imaginary pure and the function 
/ of equation (26) would become the sum of two imaginary pure terms of the 
same sign, so that it cannot be zero. In equation (30) the coefficient of the sinus 
is the sum of two non zero terms of the same sign, so that to ensure the condition 
we need to check that 

Rr)%c 



tan(coi?) = 



q + Rt] cI 



(31) 



does not holds for 770 solution of (26). The compatibility condition between (26) 
and (31) is 

_ co _ Rrjcp 
770 q + Rrjoc 2 , 

or equivalently 

(q + Rr) cl + Rr)%)co = , 
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which cannot be verified (recall that 770 7^ q, so that cq 7^ 0). The lemma is 
proved. □ 

We return to the proof of Proposition 3.1. Since the previous Lemma guar- 
antees that the Jacobian matrix J of the derivatives of F with respect to 
£ and 77 is invcrtible, we can apply the implicit function Theorem at point 
(£, 77, a) = (0,770,0). Fix £ = irjo. By Proposition 3.2 the flow defined by the 
system (27) is C 1 in the parameters, so that its derivative in a coincide with the 
first term of the Taylor expansion, denoted "J^ 1 ); 



dvI/W = 



iq -170 



dx + i 



1 

1 



(32) 



Here, denoted the solution of the deterministic problem (a = 0). Let M be 
the matrix appearing in the drift term of the above equation; the solution can 
be computed explicitly: 

*W(x)=iJ*exp(M{x-y)) ^\y)dW v , 



where 

exp(M(x-y))&°\y) 



i— cos (co(x - y)) sin(coy) 
Co 

+ 2i— sin (c (x - y)) sin(c j/) 
c o 

+ i— sin (c (x - y)) cos(c y) 
co 



i — sin(coa;) + 2i 



co 



• 



(cq(x - y)) sin(c y) 



It follows that 



9 CT F(0,7 ?0 ,0) = i 



exp 



(M(iJ-2,))*<°> 



dVF„ 



Observe that F(Q, 770, 0) is imaginary and, symmetrically, d v F(Q, 770, 0) is real. 



Set a := F{0, 770, 0) 



CO 



sin (c i?) — i?^r cos (co-R). We can write 



J- 1 = 



and from the formula 



we get that c^t? is given by (28) and cV^ = 



(0,r7o,0) 



□ 
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Corollary 3.4. When q = (2n + l)n/2R for some n G N, in the deterministic 
case we have the creation of what is sometimes called a "quiescent" soliton. 
Also in this case, the stochastic perturbation can modify, at first order, only the 
amplitude of the soliton. As a result, the "quiescent" soliton is destroyed and 
with probability 1/2 a true soliton is created. 

Proof. In the case of a "quiescent" soliton Lemma 3.3 still holds. Indeed we 
have now cq = q and qR = (2n + \)n/2, so that equation (30) becomes 



<%F (0,0,0) = — i 



Cn C 



C 

sin (cqR) + ti?-4 cos (c R) = ^fi- ^ 
c o 1 



and the determinant of the Jacobian is not zero. One has then a = ±l/g and 

8 a F(0, 0, 0) = - f sin(gi?) dW v = tWr , 
Jo 

which is real. Therefore, 

A true soliton is created whenever Wr > 0. □ 

Remark 3.5. The result of the above Remark holds in greater generality. In- 
deed, substituting the white noise with a general process Q x one obtains the 
same result: d a j] = q J Q R Q x dx, so that a true soliton is created whenever the in- 
tegral is positive. We see that for this result to hold we only need the smallncss 
assumption on the perturbing process, but this process needs not be rapidly 
oscillating. 



3.3. Small-intensity complex white noise 



The limit case of a small-amplitude complex white noise is very similar to the 
case of the real white noise treated above. Take Uq(x) = (q + aW x ) lr 0) ffl(a;) 



where W x = Wx 1 ' + iW x ^' is a complex Wiener process. We have 

Proposition 3.6. For qR > ^ and in the limit of a small complex white noise 
type stochastic perturbation of the initial condition, the velocity and amplitude 
of solitons are perturbed at first order by small, zero-mean, Gaussian random 
variables, which are given by 



-(2) 



qsm(c R) W R 2} + J R 2^f sin (c (R - y)) sin(c y)dwi 2) 



<r 
7? 



q sm(c R) W R 1] + j" 2?f sin (c Q (R - y)) sin(c 2,)d^ 



Rr] sin (cq-R) — Rj 1 cos (coi?) 



(i) 



<r 

73 



Rijo 



sin (coi?) 



R^-cos 

CO 



(coi?) 
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where cq := \/q 2 — rfo, and r]o is the amplitude of the soliton of the unperturbed 
system. 

Proof. This proof is similar to the one of Proposition 3.1. An analogous of 
Proposition 3.2 holds in this setting, and Lemma 3.3 remains unchanged (note 
that in Lemma 3.3 we work on the deterministic equation). From equation (32) 
onward one has just to remember that W is now complex. □ 

Remark 3.7. In this case the first order perturbations of the velocity and 
amplitude of the soliton have the same law and are independent. 

Remark 3.8. "Quiescent" solitons are perturbed in both amplitude and speed, 
leading to the possible creation of true solitons with nonzero velocity. The same 
result holds when perturbing the initial data with more general complex pro- 
cesses. 



4. Stability of KdV solitons 

In this section we study our second example, the KdV equation. We focus on the 
soliton components of the solution. Solitons of the KdV equation are identified 
by an imaginary number £ = it], defining both the velocity and amplitude of the 
soliton, which are related. We start again by reporting in Paragraph 4.1 some 
classical results on the background deterministic solution. In Paragraph 4.2 we 
analyze how this solution is modified by the introduction of small-amplitude 
white noise perturbation of the initial condition: the main result is contained 
in Proposition 4.1, while Proposition 4.4 deals with the case of "quiescent" 
solitons. The last paragraph deals with the case of a perturbation of the zero 
initial condition. 

4-1- Deterministic background solution 

In [Mu78], Murray obtained a C°° solution for the KdV equation with a deter- 
ministic "box-shaped" initial condition Uq = q lr-^ra (x). He showed that in 
this case the Jost coefficient a extends to an analytical function in the upper 
part of the C~plane. Only in the case of positive values of q, a has a finite number 
of zeros on the imaginary axis £ = ir\ for < r\ < yjq. 

We report some explicit computations on our similar deterministic the 
results will be used below, and study some properties of the soliton components 
of the solution. First, let us construct explicitly the solutions of the deterministic 
equation, which we call lpq. Take here U$ = q 1[o,a] {x) for some q > 0. Due to 
[AC, Lemma 2.2.2], we can assume £ = ir/; we need to solve 

J r) 2 (p x<0,x>R , . 

^ XX _ { (-q + V 2 )^ X€[0,R) ■ ^1 
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For i] — the only intcgrable solution is tp = 0. Soliton components correspond- 
ing to rj > must satisfy 

ip = Cl e ,lx x<0, (34) 
tp = c 2 e- rlx x>R. (35) 

For x £ [0, R] one can rewrite the problem as 

d tp = tp Ax 



d tp = (—q + Tj)ip dx 

Note that the sign of r] is not relevant (we take it to be positive) and that 
for rj > yjq the solution of (33) is monotone, so that it cannot be a soliton 
(which must be intcgrable). We therefore look for solutions corresponding to 
< T) < y/q. Set c = s/q — rf. Due to (34) and (35), we only need to solve (33) 
for x <G [0, R]. From (33) and the initial conditions 

(p (0) = ci , d x <po(0) = TfCi 
derived from (34), we get 

ip {x) = ae^ + Pe~^ , a =|(l-^), /?=|(l + ^), 



which is to say 



V 



ipo{x) = ci cosh(icx) — ic\ — sinh(zcx) = c\ cos(cx) H — sin(ca;) . (36) 



c 



V 



c 



We can set the global constant c\ equal to 1. Considering also the final condition 
(35), we obtain an equation for rj: 



q sin(£R) = 2£r/c 2 e- i? '' f q smjRy /q - t? 2 ) = 2c 2 T)^q - ifeT R ^ 

q cos(£i?) = (e - 77 2 )c 2e - fl " - \ g cos( J R v 4 _ 1?) = c 2 {q - 2if)e- Rr > " 

For = \Jq/2, the only possible solution is such that R\J q — rf = ir/2 + kir, 
which means that 

y/qR=(2k + l)Tr/V2 . (37) 
All other solutions can be found solving 



f( V ) := tan(i? ^fq~rf) - ^^\f = Q (38) 

q — 2rf 

for rj £ [0, ^/q)\{\/q/2\. The existence and the number of solutions for the 
above equation depend on the quantity Ry/q. Consider some fixed value of R. 
As it is shown below, for small values of q a first soliton is created with rf^ ~ 0. 
As q increases, the value of 7/ 1 ) increases too and tends to \Jq/2 as q tends 
to ir 2 /(2R 2 ). We have already found the solution for this specific value of q 
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(k = in equation (37)). For q larger than 7r 2 /(2i? 2 ), r/ 1 ) continues to grow. A 
second solution appears (r]^ = 0) when q = tt 2 /R 2 . A third solution appears 
at q = (27r/i?) 2 ; the values of 77W (corresponding to the i th soliton created) 
continuously increase as q grows, but remain ordered: rf %s> < rj^ for i > j. 
Therefore, the number of solitons created is \_R^/q/n\ + 1. 

A few examples of f(rj) are plotted below. We have taken R = 1 and different 
values of q. The first critical points (when new solitons are created) correspond 
here to q {2 ~> = tt 2 - 9,87, q^ = 4tt 2 - 39,48, q^ = 9tt 2 - 88,83, q& = 16tt 2 . 
The almost-vertical line appearing near r\ = s/q/2 for q = 44 reflects the 
fact that we are near the critical points of (37): from (37) for k — 1 we have 
q = 9tt 2 /2 - 44.41. 




Let us take a closer look at the case q 0. We assume q — q^e and look for 
the first terms of the expansion of 77 in e: r\ = 770 + r)i£ + T]2S 2 + 0(e 3 ). The first 
term of the expansion in e of the function / defined by (38) gives 

(ord. 0) /(t?) = t a n(R^~n5) + + 0(e). 

Since the order-zero term cannot be made equal to zero, we must skip the first 
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term in the expansion of r\: r\ = t]\e + 0(e 2 ). Looking at equation (38) at first 
order 

(ord. 1) f(ji) =(Rqo- 2m) + 0(e 3 / 2 ) = 

V / Via 

we obtain rji = -^j 2 -. Pushing the expansion further, one can obtain the following 
order coefficients. At order two we have: 

(ord. 2) /fa) = - (±i? 3 ql + 2 ??2 ) tl + (e 5 / 2 ). 

One has then 

V=^e-^ + 0(s% (39) 
showing that in the limit n is of the same order of q. 



4-2. Small-amplitude random perturbation with q > 

We want now to add some random perturbation to this "box-shaped" initial 
profile. 

Let Uq = (q + <jW x ) lroffl(x) be the initial condition of the KdV equation, 
where W is a standard Wiener process. Since solitons correspond to zeros of the 
complex extension of a, which in turn must be located on the imaginary axis, 
we look for bounded solutions of equation (7) for Q = it], r\ € R. The first main 
result is contained in the following Proposition. 

Proposition 4.1. In the limit of a small white noise type stochastic perturbation 
of the initial data, the parameter r) defining the velocity and amplitude of the 
generated soliton is perturbed at first order by a small, zero mean, Gaussian 
random variable, which is given by 



/,f <p {R - x)ip (x)dW x 



cos(co-R) 



rjoR - 



sin(coi?) 



3)X>+2?)g-R 

Co 



'/,; 

7HT 



(40) 



Here, (p is the deterministic solution of (33), given by (36), cq := \J q 2 — rj 2 , 
and t]q is the amplitude of the soliton of the unperturbed system. 

Proof. We need to solve 



dip = ipdx 

dip = {—q + i] 2 )tp dx + aip o dW x 



(41) 



for x £ [0,R] with initial conditions ip(0) = 1, <J5(0) = 77. An analog of Propo- 
sition 3.2 holds, and it provides the uniqueness and regularity results for the 
solution. 
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Proposition 4.2. The stochastic differential equation (4-1) defines a stochastic 
flow ^x'^ = ((p x ,tp x ) T of C 1 -diffeomorphisms, which is C 1 also in the param- 
eters r),o~. 

As we did for the NLS equation, we look at the flow provided by the above 
Proposition 4.2 and define a function of the flow at the point x = R as F(rj, a) := 
(f(R) + r]ip(R): this remains a function of the two parameters (77,(7). To prove 
Proposition 4.1 we use the implicit function Theorem to show that F(r],o~) has 
a unique zero in some open set containing (770,0), the point corresponding to 
the deterministic solution. The following Lemma ensures that the hypothesis of 
the implicit function Theorem are satisfied. 

Lemma 4.3. Let F(j), cr) be the function defined above. Then, for all 770 corre- 
sponding to soliton components of solutions of the deterministic problem, d v F(r]o,0) ^ 
0. 



Proof of Lemma 4-3. For 77 
that c = i/g-??o) 



r)o and a = we have that (p = ipo and (recall 



d v <p (R) = sin(co-R) 



l + 77 i? , Vo] RVo 



— Y- cos(cq-R) , 



dr,p(R) = cos(c i?) [1 + #770] +^-sm(c R)[l + T] R] = [l + rj R]cp (R) 



Since 



d n F 



d v (p + ip + rjd v (p , 



d v F(r] O ,0) = cos(cq-R) 



VoR 2~ 

4 



sin(co-R) 



3?7o 



4 



CO 



The coefficient of the sinus is strictly positive (the coefficient of the cosinus has 
instead at least one zero for 770 G [0,y^], since it is positive for 770 = and 
negative for 770 — > yfq)- Therefore, we only need to verify that the following 
equation 



g{R,q,Vo) ■= tan(c i?) 



2 + VoR 



co 



< o 



(42) 



is not satisfied, knowing that 770 is a value corresponding to a soliton solution 
of the deterministic equation. As we have seen, we can either have 770 = \J q/2 
if condition (37) is satisfied, or else 770 is given as the solution of equation (38) 
in (0,^)\{v^72}. 

In the first case we have that co = \/<z/2, so that condition (37) implies that 
cos(coi?) = and sin(co-R) = ±1. Therefore, d v F 7^ 0. 
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Consider now the second case. We look for points 77 G (0, ^/q)\{\/q/2} such that 
fiv) = 9(v) = 0- If such a point exists, then 



which also reads 



or equivalently 



2-qc 
q — 2r\ 2 

2 V 



q-2r) 2 



2 + r,R-*£ 

r. 



(2 + 7?fl)c 2 - ifR 
(377 + 2r) 2 R)c 2 + T] 3 



= 



{2-q 2 + 2r; 3 fl)(g - if) + 2t? 4 (1 + r)R) + (2 + 7?fl)(g - t? 2 )^ - ifRq 
(q - 2t 1 2 ){(3i 1 + 2rj 2 R){q - r, 2 ) + r? 3 ] 

For rj < s/q/2 the denominator is positive and the numerator 

(2?7 2 + 2ifR)(q - rj 2 ) + 2t7 4 (1 + r]R) + (2 + V R){q - V 2 )q - rfRq 
> r)R(q - if) - ifRq > , 



(43) 



so that the fraction (43) is positive and cannot be zero. For ^Jq/2 < rj < ^fq 
the denominator is negative and the numerator is larger than 2rf J R — rj 3 Rq = 
ifR{if — q) > 0, so that the fraction (43) is negative and cannot be zero. The 
Lemma is proved. □ 

We return to the proof of Proposition 4.1. Since the above lemma guarantees 
that d v F(rjo, 0) ^ 0, we can apply the implicit function Theorem. We have (at 

m) 

dd a (p — d a (p dx 

dd a fi = (-q + r] 2 )d a (pdx - (ip + <jd a (p) o dW x 

By Proposition 4.2 the flow $i ?) ' <T ' > is C 1 in the parameters, so that its derivative 
in a at a = coincide with the first term of the Taylor expansion, which is the 
solution of 

dd a ip = d a (pdx 

dd a (p = {-q + r]^)d a ipdx - ip dW x . 



In matricial notation 
d 

with 

M = 



d a <P 



d$ 



-c 2 



$dx 



( ° 







z 2 



Mx 



cos(cx) - sin(ca;) 
-c sin(ca;) cos(cx) 
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The solution is 



= $(0) - f 
Jo 



,M(x-y) 



We have 



d„<p(x) = 
d a <p(x) = 



sin (c(x — y)) cos(cy) H — sin(cy) 

v ' L c 

os (c(x — y)) cos(cy) H — sin(cy) 



dW„ 



dW, 



1/ • 



Therefore, 



9 ct F(t7o, 0) = 9 CT ^(i?) + VodM^) = - / MR ~ a?)^)(x)dW x 

Jo 

and we obtain, at first order, the perturbation of the value of the parameter r\ 
defining the velocity and amplitude of the soliton: for every (q,R,rjo), 



d a n(a) 



d a F 



Jo 


cos(c (R-x)) + ^sin(co(ii-a;)) 


cos(coa;)+ ^ a sin(cpa;) 


dW x 


cos(co-R) 


2 + m R - 

c o _ 


+ sin(coi?) 


c ~cj 





The Proposition is proved. 



□ 



We turn now to consider the case when the stochastic perturbation can result 
in the creation of a new soliton. As for the NLS equation, this happens only 
for specific "critical" values of q (and R). The result is presented in the next 
Proposition. 

Proposition 4.4. If we are at a critical point, which is to say yfqR = nir for 
n£No, a small amplitude white noise type stochastic perturbation of the poten- 
tial may create a new small-amplitude soliton. The condition for the creation of 
a new soliton is that the zero-mean, Gaussian random variable (4-5) is positive. 



Proof. The deterministic background solution is (c = yg — rj 2 ) 



ipo = cos(cx) H — sin(cx) 
c 



(44) 



We have Uq = q + uW x ; we want to apply the implicit function Theorem to the 
function F(r],a) defined above, at the point (0,0). At this point 

d v ip{R) = -L S m(^qR)=0 , 
d v &(R) =cos(y/q~R) = ±1 . 
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Therefore 

d v F[0,0) = d v tp + tp + r)d n <p = 2cos(Vgi?) = ±2 , 
and the implicit function Theorem is applicable. We also have (again at r\ = 0) 
dda-cp — d a tf dx 

dd a (p = —q da-ip dx — (ip + ad a ip) o dW x 

As seen above, the solution of the above system coincide at <x = with the 
solution of 

{dd a ip = d a ip dx 
dd a (p = -q d a ip dx - p dW x ' 

which is given by (here, c = y/q) 

1 f x 

d a (p(x) = / sin [y/q(x - y)) cos{y/qy)dW y , 

\f1 Jo 



d a ip{x) = cos (y/q{x - y)) cos(y/qy)dW y . 

Jo 

We have therefore 

r R 

d„F(0,0) = - / coa(y/q(R-x))coa(y/qx)dW x , 
Jo 

and since cos(y/qR) ± 1 

fl V /nn , J R co S (^q(R-y))cos(^qy)dW v 



2 



n 



d v F y ' ' 2cos(^/gi?) 
cos 2 (^x)d^ . (45) 



In this case, a new small-amplitude soliton, corresponding to i] = av, is created 
whenever v > 0. □ 



4-3. Small-amplitude random perturbations with q—0 

We analyze now the case in which the initial condition is the pure stochastic 
perturbation: contrarily to the NLS equation, even this small initial condition 
can generate a soliton. In this setting there is no (non trivial) solution in the 
deterministic case. This case is obtained at the critical point y/qR = 0, where 
the first quiescent "soliton" is created. 

Proposition 4.5. For q — ; a small amplitude stochastic perturbation of the 
potential may create a new small- amplitude soliton. For a white noise type 
perturbation W x , the probability that a new soliton is created is 1/2 and the 
created soliton corresponds to the value rj = aWfj. 
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Proof. Again, we want to use the implicit function Theorem; take the limit of 
the deterministic solution (44) as q — > 0: 



We have 



and (rj = 0) 



Therefore, 



ip = 1 , <p = 



d n F(0, 0) = <9, ; <^ + 99 + Tjd T i^p = 1 



d</? = t/j da; 

d<^ = — a tp o dW x 



dd a <p = do-ipdx 

dd a (p = -(tp + od a ip) o dW x 



d a F(0, 0) = d a lp + r]d a ip = - [ ipdW x 

Jo 



-Wr 



It follows that 



and a single soliton corresponding to r\ = oWr is generated whenever Wr > 0, 
which means with probability |. □ 

Remark 4.6. The same result of Proposition 4.5 holds with more general pro- 
cesses: if, instead of a white noise, we take Uo(x) = a Q x in (33), where Q x is a 
generic stochastic process, we have (at rj = 0) 

dd a tp = d a (p dx 
dd a tp = —Q x Up + ad a ip)dx ' 



d a F(Q, 0) = d a ip + r\d a ip = - I Q x dx 

Jo 



and 

cR 



In this case, a single soliton is created whenever the noise introduced has positive 
mean, and it corresponds to rj = a J Q R Q x dx. 
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Remark 4.7. The mass of a soliton of the KdV equation is M v = J R U{x) dx = 
477. Here, we have introduced a perturbation of mass Mjj = a J R Q x dx = 77. 
We see therefore that the soliton created has a larger mass that the initial 
perturbation, implying that the radiative part (going in the direction opposite 
to the one of the soliton) has absorbed a total mass of 377. 

The energy conversion efficiency from a small noise source to a soliton is 
very poor, since the input energy is of order a 2 (Ejj = a 2 J R Q 2 dx) , while 
the energy of the created soliton is only E v = J R U 2 (x) dx = ^-rf ~ a 3 , for a 
smooth source. For a white noise source, the input energy is infinite, while the 
energy of the created soliton is finite. 



5. Conclusion 

For the NLS and KdV equations, the study of soliton emergence from a local- 
ized, bounded initial condition perturbed by a wide class of rapidly oscillating 
random processes can be reduced to the study of a canonical system of SDEs, 
formally corresponding to the white noise perturbation of the initial condition. 
The integrated covariance is the only parameter of the perturbation process 
that influences the limit system of SDEs. From the study of this limit system, 
one obtains quantitative information on the modification of solitons due to the 
random perturbation. 

For the NLS equation we have a threshold effect: if for the deterministic initial 
condition the integral L Uo{x)dx exceeds n/2 at least one soliton is created. In 
this case, a small amplitude random perturbation of the initial condition results 
in a small variation in amplitude for the created soliton. However, the speed 
of the created soliton is not substantially modified, unless the perturbation is 
complex. 

On the other hand, since the KdV equation does not present a threshold 
phenomenon, a small amplitude random perturbation always results in a small 
variation in both speed and amplitude of the created soliton, for any non nega- 
tive initial condition. 

As for "quiescent" soliton, for both NLS and KdV a stochastic perturba- 
tion has a positive probability (depending on the type of perturbation used) of 
creating a new real soliton. 
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